Local minimizers and low energy paths in a model 
of material microstructure with a surface energy 

term* 



A family of integral functionals T which, in a simplified way, model 
material microstructure occupying a two-dimensional domain f2 and 
which take account of surface energy and a variable well-depth is 
studied. It is shown that there is a critical well-depth, whose scal- 
ing with the surface energy density and domain dimensions is given, 
below which the state w = is the global minimizer of a typical F 
in T . It is also shown that m = is a strict local minimizer of 
F in the sense that if u 7^ is admissible and either llwHiSfQ-) or 
S : \vy\[x^y) > 1}) is sufficiently small (with quantita- 
tive bounds given in terms of the parameters appearing in the energy 
functional F) then F{v) > F{0). Low energy paths between u = and 
the global minimizer (in the case of a sufficiently large well-depth) are 
given such that the cost of introducing sets {{x, y) £ V, : \vy{x, y)\ > 1} 
of positive measure into the domain Q may be made arbitrarily small. 

1 Introduction 

The energy functionals we shall consider in this paper are related to the one 
used in jKM94 j but with some important differences. The original Kohn- 
Miiller functional is 
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where Q = [0, L] x [0, 1], L > 0, e > is a small parameter which is sometimes 
referred to as the surface energy density, and dx is shorthand for dC'^{x). 

Provided suitable boundary conditions are imposed, the global mini- 
mizers of Ekm model in a simplified way the fine scale microstructures 
that are observed to some degree at austenite-martensite interfaces in shape 
memory alloys. The second derivative term tempers the oscillations in the 
y— direction that any globally minimizing sequence will develop. Using a 
subtle argument Kohn and Miiller showed that the global minimization can 
be viewed as a straightforward competition between the term in ti^, effec- 
tively a measure of 'spread', and a version of surface energy derived from 
the terms in Uyy dx and f^iuy — 1)^ dx, among functions u with \uy \ = 1 
a.e. The result is that the infimum of the energy scales in e as though it 
were evaluated at the now well-known branched microstructure (see [KM94J 
for details). 

However, the functional Ekm is less useful in understanding the role of 
= as a local minimizer and other details of the energy landscape. To 
this end we introduce new functionals Ei{-] e, A), E2{-; e. A) and E^i^-; e. A) 
based on Ekm but with some extra features. The family of functionals 
mentioned above consists of £i(-;e,A), E2{-;e,A) and -E3(-;e, A) defined 
below in (jl.ip . ()1.2p and (II. 3p respectively as the parameters e and A vary. 

Definition 1.1. 

^i(n;e,A) = / e'^ul + uldx + AC'^{A{u)) (1.1) 
Jn 

E2iu;e,A) = [ e^\Viuy)\'^ + uldx + AC\A{u)) (1.2) 
Jn 

E3{u;e,A) = [ e^\V{V{u))\^ + uldx + AC\A{u)) (1.3) 
Jn 

B{u) = {(x,y) G : \uy{x,y)\ > 1} 
A{u) = n\B{u). 

In each case the functionals take the form 

Ei{u) = Si{u)+ [ WA(Vu)dx, 
Jn 

where 

WA{a,b) = + Ax^-i,i){b) 

and where Si{u) is a surface energy term. In particular, S3{u) will henceforth 
be written in the more conventional way 

Ssiu) = [ |D\pdx. 
Jn 
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In the following we suppress the dependence of the Ei on e and A for 
brevity. We study the behaviour of each Ei in the class Ai of admissible 
maps defined below. First we define the subclass of l^p'^(^7,M) C W^''^{n,R) 
of functions satisfying periodic boundary conditions (in the sense of trace) 
at the top and bottom of the domain by 

W^'^{n, R) = {ue W^'^{n, M) : u{x, 1) = u{x, 0) < a; < 1}. (1.4) 

Then 

Ai = {uGW^'^{n,R) : Uyy e L^{n;R),u{0,y) =00<y <l} (1.5) 
A2 = {u G Wp^'2(0,M) : G P^^'2(O;R),u(0,y) = < 2/ < 1}(1.6) 
^3 = {ueW^^^{n,R):VueW^''^{n;R),u{0,y)=00<y<lll.7) 

The new features of these models relative to the original Kohn-Miiller 
functional are summarized here and discussed below. They are: 

(i) a variable well-depth A (see below for its definition); 

(ii) a convex potential in a neighbourhood of (0,0); 

(iii) the possibility of a cost, which may be zero, associated with the ap- 
pearance of sets B{u) of positive measure in Q, analagous to a lower 
bound on the cost of 'nucleation' of martensite in austenite. 

The term C'^{A{u)) is intended to mimic the behaviour of the term 
J^{Uy — l)^dx appearing in the Kohn-Miiller functional Ekm in the fol- 
lowing sense. In order for Ei{u) to approach its infimum \uy\ < 1 can only 
occur on a set of small measure: 'most' values of \uy\ will be near to or 
larger than 1. Looking at Ekm and referring to the argument they give we 
see that most values of the gradient Uy of the global minimizer in that case 
will be near ±1. The price we pay for replacing J^iuy — l)^dx with a term 
proportional to £?{A{u)) is that large values of \uy\ are not penalized as they 
would be in Ekm- But this turns out not to matter a great deal, as will 
become clear later. In fact, the nature of the global minimizer of Ei in Ai 
in a scaling sense can be deduced from the Kohn-Miiller argument when A 
is large enough, although we do not pursue this in the present work. When 
A is in the range (0, CeL~^) for an appropriate dimensionless constant C 
it happens that u = is the global minimizer of Ei in Ai for i = 1,2,3. 
This behaviour with respect to varying the well-depth A may simply be an 
artefact of the low dimensions and the choice of boundary conditions. In 
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any case it will be necessary to know just how large A needs to be before 
u = ceases to be the global minimizer of Ei in Ai for i = 1,2, 3. 
The idea of introducing a well-depth 

A = Wa(0,0)-Wa(0,±1) 

comes from the Ball- James theory of martensitic phase transformations. The 
theory asserts that the stored-energy potential should change in a certain 
way as the temperature changes; we synthesize this by varying A, with A = 
corresponding to high temperature stored energy function and A = 1, say, 
to low temperature. The global minimizer in the case A = is n = which, 
in the full three-dimensional models, would be referred to as Austenite. See 
|BJ87| . |BJ92j for further details. We do not attempt to introduce dynamics. 
It is shown in Section[3]that ti = is always a local minimizer of Ei regardless 
of the size of A; the only effect A has on local minimality is through the 
size of the neighbourhood A/i, say, of u = in Ai on which u = satisfies 
Ei{v) > Ei{0) for all v in Afi for i = 1,2, 3. See Section^ and in particular 
Theorem l3.H for details. When the well-depth is large enough it also appears 
in the scaling of the energy associated with the global minimizer; we do not 
pusue this in the present work. We remark that in the case of Ei one could 
follow the method given in |KM94j with only minor changes. 
Note that the potential 

WA{a,b) = + Axi-i,i){b) 

is convex in a neighbourhood of (a, b) = (0, 0). This is sufficient to establish 
that ti = is an L^-local minimizer of Ei in Ai- It is not necessary, though, 
as examples of Taheri show |Ta021 Section 4]. In fact the potential Wa 
can be bounded below by a strongly convex potential in a neighbourhood of 
zero by 'borrowing' some surface energy and applying a suitable Poincare 
inequality. For example, for any A € (0, 1), 

Ea{u) > / (1 - A)e2|Vny|2 + CnXe^l + t^A(Vu) dx; 
Jn 

so that WA{s,t) := WA{s,t) + C^Ae^t^ is strongly convex in a neighbour- 
hood of (0, 0). If we are allowed to vary the domain height, for example by 
taking 

Qh = [0, L] X [0, h] 

one can use the same procedure to bound the Kohn-Miiller potential below 
by a potential that is strongly convex at (0, 0). The reason this works is that 
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the Poincare constant Cq,^, say, becomes large when h is small. See Taheri 
|Ta02j for other interesting examples. When h = 1 the Poincare constant 
is not large enough for this trick to work, from which it follows easily that 
= is not a local minimizer of Ek m in ^ (provided e is sufficiently small) . 
We note that the idea of combining the strong convexity of W with higher 
order terms in order to guarantee local minimality has been studied before 
[Ta02j . |Ba05j . Although we don't use these arguments directly they are, in 
view of the comments above, one of the main reasons that we can expect 
li = to be a local minimizer of Ei in Ai for i = 1,2,3. 

The results of Section [3] are based on an apparently new inequality which 
relates all three terms appearing in Ei and which, together with more stan- 
dard estimates, yields explicit functions r(e,A) and s(e,A) (given in (jS.ip 
and (13. 2p respectively) such that 

\H\L^(n) < r{e,A) or C^Biu)) < s{e,A) =^ E^{u) > ^a(O), (1.8) 

provided u 0. It is doubtful whether the scalings involved are optimal 
for reasons explained in Section [3l Even so, it is still a stronger and more 
explicit result than L^-local minimality. The results for the functional Ei 
are easily carried over to the functionals E2 and E^. 

We discuss in Section [J] the effect of modifying the surface energy term. 
The reasoning set out in Section [4] points out that any path that connects 
u = with the global minimizer must, provided A is large enough, pass 
through a state at which sets B{u) of positive measure first appear. This 
is the basis for a calculation which tests whether such states automatically 
cost a certain minimum amount of energy, analagous to a 'nucleation cost'. 
The answer is that sets B of positive measure can appear at arbitrarily small 
energy (measured by any of the Ei). It may help to recall that for a given 
function u in Ai the set B{u) = {{x,y) € Q : \vy{x,y)\ > 1}. We surmise 
that if there is to be a nucleation cost then it must be as a consequence of 
some factor beyond the details of the energy functional itself, such as the 
dynamics, for example. See SectionS]for further details. The sets B{v) also 
play an important role in the local mimimality results of Section [3l This is 
another reason that the introduction of sets B{u) of positive measure is of 
interest to us. 

1.1 Notation 

One and two dimensional Lebesgue measure are denoted respectively by 
and C"^ throughout the paper. The usual Sobolev space notation W^''^ is 
used to denote the class of A;— times weakly differentiable functions all 
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of whose derivatives lie in L^{^}). The norm of a function v is denoted 
by \\v\\p or depending on the context. The class of absolutely con- 

tinuous functions is denoted AC, and those functions which are absolutely 
continuous along almost all lines parallel to the coordinate axes are written 
ACL. See [Zi] for further details. The characteristic function of any set S 
is written xs- All other notation is standard, with the possible exception 
that the value of the dimensionless positive constants C appearing in var- 
ious inequalities may, where no confusion arises, change from line to line. 
Where it is necessary to distinguish between positive constants we shall use 
the convention that if C and c appear in the same calculation then c < C. 
We will also employ the convention that roman letters x represent vectors 
in M2. 



2 The effect of a variable well-depth 

When A = it is clear by inspection that n = is a global minimizer of Ei 
in Ai for i = 1,2,3. Therefore the following definition makes sense: 

Definition 2.1. For each i = 1,2,3, 

A, = sup{A > : u = globally minimizes Ei in Ai}. 

Each Ai will be referred to as a critical well-depth. The ordering 

Ei{v) < E2{v) < Esiv) (2.1) 

for all appropriate v, together with the fact that -E'i(O) = -£2(0) = -^3(0) for 
fixed e and A, implies the inclusions 

{A : Ei{v) > Ei{0) yv e Ai} C {A : E2{v) > E2{0) \fv £ A2} 
{A : E2{v) > ^2(0) yv e A2} C {A : E^iv) > E^iO) Mv € ^3}, 

and hence 

Ai < A2 < A3. (2.2) 

It will be shown in this section that all the Aj scale alike with respect to e 
and L in the sense that there are dimensionless constants c < C such that 

ce . . . Ce 
- < Ai < A2 < A3 < — 

provided e is sufficiently small. 



Local minimizers and low energy paths 



7 



To begin with, Proposition 12.21 below shows that there is a constant C 
depending only L such that Ai > Ce^. This follows relatively simply by 
using only the surface energy term in the functional Ei and it turns out 
to be a crude lower bound on Ai. A better (i.e., larger) lower bound is 
obtained in Section 12.21 where it shown that 



The optimality (in terms of its scaling in e and L) of this lower bound is 
proved by evaluating the energy of a particular element of ^3 in sec- 
tion 12. 3[ The structure of is given in section 12.31 It is not a branching 
microstructure, unlike the global minimizer obtained when Aj ~ 1 and cor- 
responding to the model studied by Kohn and Miiller. The main result of 
this part of the paper is: 

Theorem 2.1. With Ei, Ai and Aj as above, there are dimensionless con- 
stants c < C and eo > such that 



Moreover, for all A the global minimizer Ui of Ei in Ai exists, and when 
is large enough it satisfies C'^{B{Ui)) > 0. 

The proof of Theorem 12.11 rests on Propositions 12.31 and 12.51 in Sections 
12.21 and 12.31 below. The lower bound is proved with the help of an interpo- 
lation inequality; the upper bound corresponds to the statement concerning 
the optimality (in a scaling sense) of the lower bound. The proof of the 
existence of the global minimizer Ui of Ei in Ai is a relatively straightfor- 
ward application of the direct method of the calculus of variations. It is 
given in an appendix for completeness. The last assertion of Theorem 12.11 
can be deduced from the comparison function constructed during the proof 
of Proposition 12.51 

2.1 A simple lower bound on Ai 

We recall that 



A(u) = {{x,y) G n : \uy{x,y)\ < 1}. 

The following lemma will be used to show that Ai is strictly positive; it 
will also play an important role in Section [3l 



^ < Ai < — 2/ < e < eo- 




where 
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Lemma 2.1. Let u £ Ai satisfy C'^{B{u)) > and define for each x G [0, L] 

I, = {x}x[0,l] (2.3) 
U{B) = {x e[0,L]: C^{l^nB{u)) >0}. (2.4) 



Then 



In '^yy 



where 



£2(5(u)) - r(l-r) 

■ ^Hn(i3(n)))- 



> Z7^^ (2-5) 



Proof. We begin by remarking that the quantity on the right-hand side of 
(j2.5p is finite under the assumptions of the lemma. Let B = B{u) for brevity. 
Clearly C'^{B) ^ implies r > 0. By Fubini's Theorem we have that r < 1; 
the membership of u in Ai further implies r < 1. If this were not so then 
for almost all x G [0, L] the lines Ix would satisfy 

C^{lxr\B) = 1. 



Then, since Uy is absolutely continuous along almost all lines l^, we must 
have for almost all x that either 



Uy{x,y) > 1 for a.e. y € [0,1] 



or 

Uy{x,y) < —1 for a.e. y G [0, 1], 

contradicting the periodic boundary conditions imposed along y = and 
y = I. In fact, this argument shows that C^{lx Ci B) < 1 for almost all x. 

Now we proceed with the proof of inequality (j2.5p . Let x G II(i?). The 
argument above shows that we may assume there is at least one open set 
3^1 C Ix^B on which Uy{x,-) > 1 and at least one other open set 3^2 C Ixf^B 
on which Uy < —1. We may suppose that yi := sup3^i < inf3^2 =: 2/2, so 
that the intervening set is [yi , 1/2] • It is easy to check that the minimum of 
the scalar functional 

ry2 

f - / iff dy 

Jyi 

among / G W'^^'^{[yi,y2],^) satisfying f'{yi) > 1 and /'(y2) < -1 is ^^4^- 
The minimization calculation is of the 'free endpoint' kind, so that it pre- 
scribes optimal values for the differences f{y*) — f{yi) and f{y2) — f{y*) 
where yi < y* < 2/2 satisfies f'{y*) = 0. Note that y* exists because 
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y (yi+y2) \'^ 



Vy S AC (/a;) by assumption. The minimizer is f{y) = - — y^^y^ i with 
y* = Since this calculation is elementary we omit the proof. 

The result is 

Uyy\x,y)dy > 



y2 - yi 



But since 



we must have 



y2-yi < C\AnQ 



Integrating over x G n(i?) and applying Jensen's inequality gives 
u„„ dx > 



Dividing by C^{B) and rearranging yields inequality (|2.5p . □ 

Remark 2.1. The global minimizer v of the functional on the left-hand side 
of (|2.5p is such that 7" = |. Though easy to construct, v can never belong 
to Ai because it violates the boundary condition n = at x = 0. 

In the rest of the paper it will be useful to have a label for those elements 
u of Ai for which C'^{B{u)) > holds. 

Definition 2.2. Af := {u e Ai : C^{B{u)) > 0}. 

Proposition 2.2. Ai > C2e^ > Cie^ > 0, where 

Proof. By Chebychev's inequality, C2 > Ci. Therefore we need only prove 
^1 > C2e^. By Remark 1 2 . 1 1 ab ove C2 = 16, but the infimum is not attained. 
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Now 

> {16e^ - A) C\B). 

So if A < 16e^ then « = is a global minimizer of Ei, and hence Ai > 
16e2. □ 

2.2 A refined lower bound on Ai 

In this section we show that there is a dimensionless constant c such that 
Ai > ^ for all e. This improves on (i.e. increases) the lower bound obtained 
in Section [2.11 The reason for the improvement is essentially that the term 
dx is brought into play. 
We shall make use of the standard interpolation inequality 

\ C fyy dy + [\f - pf dy>C C dy, (2.8) 
^ Jo Jo Jo 

which holds for some C > 0, all p, all non-zero a and all / € VF'^'^([0, 1], M). 
(See, e.g., [GTl Section 7.12].) Let u e Ai, x e n{B{u)) and take p = 0, 
f{y) = u{x,y) in (j2.8p above. Using the inequality 

1 

ulix,y)dy>C\l,nB{u)), 

integrating over x € n(i?) and using Fubini's Theorem we obtain 

^ / ul dx + a^ [ u2dx>C£2(s(n)) (2.9) 
^ Jn JQ 

Minimizing the left-hand side of (|2.9p over non-zero a we see that 



^4,clxJ ^J^u',.) >™z^. (..10) 

Note that the constant C is independent of the dimensions of the domain 

n. 



( l\.2.\KcC\B{u)) 
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We also need the standard Poincare inequality 

C f 

uldx>-^ / n^dx, (2.11) 

which uses the boundary condition u = along x = 0. The constant C is 
independent of the domain dimensions. 

Proposition 2.3. There is a dimensionless constant c > such that Ai > 
^ for all e > 0. In particular, the lower bound on Ai stated in Theorem \2.1\ 
holds. 

Proof. Let u ^ Ai and set B = B (n) . By definition of Ei and from inequal- 
ities (fZTO]) and (plT]) we have 



Letting 



/n^'dx 



C^^{B) 

we see that the right-hand side of (I2.12P above has the form 

'Ce^ Ct 



4t 

the term in brackets is minimized when t = ceL for some constant c. Prom 
this it follows that any A < ^ is such that li = is a global minimizer of 
Ea. Therefore Ai > f . □ 

In some cases one can do better than Proposition 12.31 The following 
lemma shows that the lower bound on Ai obtained above is correct with 
constant c = 1 provided condition (12.131) below holds. This supplementary 
condition amounts to a strengthening of the boundary condition along y = 
and y = 1; it is satisfied, for example, by all admissible functions having 
compact support in Q. 

Lemma 2.2. Let v G Af satisfy 

{VyV:c)ydy = (2.13) 
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for almost every x in [0, L]. Then for almost every x in n(I?) 

/ €^vly + vld:^>eC^{l^f\B). 
Jn 



(2.14) 



From this it follows that Ai > . 
Furthermore, if v inA2 satisfies 



A 



C\B) 



(2.15) 



then 



E2{v) > ^2(0). 
Proof. First fix x € n(i?) for which (j2.13p holds. Then 



e^vly + w^dx > 



[0,a:]x[0,l] 



2 2 I 2 T 

e Vyy + dx 



> 2e 
= 2e 

= 2e 



VyyWx dx 



{VyVx)y - VyVyx dX 



{vyVx)ydy 



] dy \ dx 



vldy 



> e£\l^nB), 

where we have apphed (j2.13p and the boundary condition v{0,y) = for 
< y < 1 to pass from the fourth to the fifth hne. Integrating both 
sides of inequaUty (j2.14p over Il{B), dividing by C^{Il{B)) and inserting the 
resulting expression into the definition of Ei{v) gives 



Eiiv)-Ei{0) > 



£i(n(i?)) 



A C'{B), 



from which the inequality Ai > follows easily. 

Inserting the integrated from of (|2.14p into E2{v) gives 

E2{v) - E2{0) >e^J^ vl, dx + (;^Tg^ " ^) ^'(^)- 
Therefore ()2.15p implies E2{v) > E2{0) as required. 



□ 
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Remark 2.4. Any A satisfying A < forces (j2.15p to hold. Therefore 
inequality (j2.15p provides a short-cut to the proof that A2 > |; whenever 
(lOk is true. 

2.3 A sharp upper bound on A3 

We show in this section that there is a constant C independent of 0, and e 
such that A3 < ^ if e is sufficiently small. The idea of the proof can be 
explained as follows. Let us suppose that for each e > there is an element 
of A3 with the properties that 

/ e2|Z)V|2 + (O^dx < Cie (2.16) 
Jn 

C\B{v')) > C2L. (2.17) 

The constants Ci and C2 should not depend on e or L. Let A < A3. Then 
in particular 

E^iv') - ^3(0) > 
on the one hand; and, using (j2.16|) and (|2.17|) above, 

Cie- AC2L > E3{v')-Es{0) 

on the other. Thus A < whenever A < A3. Letting A — > A3 yields the 
desired upper bound. It remains to prove the existence of a map G A3 
with the properties (j2.16p and (j2.17p . 

Proposition 2.5. There exists a map in the class A3 and dimensionless 
constants Ci and C2 such that (j2.16p and ()2.17p hold. In particular, there 
is a dimensionless constant C such that 

^ - L 

whenever e is sufficiently small, proving the upper hound on A3 stated in 
Theorem \2.1\ 

Proof. Let k,l,h>0 and define : [0, Z] ^ M by 

H{x) = h — kx 
where kl = ^. Define the function w on [0, /] x [0, 2h] by 
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w{x,y) 



2H(x) 

H(x) 

y — — 



2h 



2 
2 



2H(x) 



ifO<x</, 0<y< H{x) 

if < X < /, H{x) <y <2h- H{x) 

if < X < 2h- H{x) <y <2h. 



Now extend w to [0, /] x [0, 4/i] by reflection in the line y = 2h, namely 
y) = w{x, 4:h — y) if < x < I, 2h < y < Ah. 



It can be checked that 

r4:h fl 







D^u;^ dx 





wl dx 



(ci + C2/C^ + Cs/C^) 



c^k Ih, 



(2.18) 
(2.19) 



where the q are positive dimensionless constants whose precise values are 
not important. Suppose h is chosen so that := ^ is a positive integer. 
Extend w by periodicity to [0, /] x [0, 1] and label the resulting function w 
again. A computation using (j2.18p above together with kl = !^ gives: 



1 A 



^0 



1 



D^wl^dx = 7^ C1 + C2 - +C3 - 



1 A 



Jo 



w1 dx 



C4- 



I 



(2.20) 



(2.211 



Clearly w is not an element of A because it doesn't satisfy the boundary 
condition at x = 0. But we can interpolate between w(0, y) and the function 
y 1-^ as follows. Define v : [0, 1] x [0, 1] ^ M by v{x, y) = j'w{0, y) and 
compute directly. 



1 fl /j2 

vldx = C6 — . (2.23) 



Jo 



By construction v and w depend only on the parameters h and /. The 
last step is to glue v and w together to give an element v'^ of A. Define 
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v^{x,y) 



v{x,y;h,l) 0<x<l, 0<y<l 
w{x -l,y;h,l) I < x < 21, < y < 1. 

It is straightforward to check that 

C\B{vn) = ^. (2.24) 

and that 

Finally we use (p:20]) - (p:2i]) to compute 




E3{v')-E3{0) = T^lci + caly) +C3(yl + C4/ + C5 ^ 



The domain of f"^ is Q, provided we choose I = ^. Choosing — = c e, 
where c is such that (4c(eL)2)~^ g N and |c — 1[ is minimized, and inserting 
into the above gives 



e / e / e 



E3{v')-E3{0) = ^ (^ci+C2-^ + C3(^-jVc4L + C5ej (2.25) 

3 1 AL , , 

+ cee-2L-2 +C7e - —. (2.26) 
lb 

By ignoring the term in A it can immediately be seen that (j2.16p is satisfied. 
From (j2.24p we have C'^{B{v^)) = ^, so (12.17j) holds. By the reasoning set 
out in the lines following (|2.16|) and (|2.17|) this concludes the proof. □ 

Remark 2.6. The proof of the upper bound on A3 can be obtained directly 
from (j2.25p as follows. Simply note that the inequality £'3(1;'^) — £"3(0) > 
holds because A < A3 has been assumed. In view of (j2.25p this gives 
A < ^. Hence A3 < ^. 



3 Austenite as a local minimizer of Ei 

We saw earlier that u = is a global minimizer of Ei in Ai provided A < A j . 
In this section we apply the interpolation inequality (j2.8p to the case A > Aj 
with the aim of proving Theorem 13.11 below. In fact, it suffices to prove 
the theorem for the functional Ei in view of the ordering (j2.ip and since 
£1(0) = £2(0) = -£3(0). The result is then automatically true for the other 
Ei. 
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linear interpolation 




L 



Figure 1: The construction of the test function v''. |fy| = 1 in the shaded 
regions; \Vy\ < 1 in all other areas. The basic repeating unit shown with a 

bold outline in the region ■§ < x < L has a vertical lengthscale of order . 



Theorem 3.1. For each e > and A > it is the case that Ei{v) > Ei(0) 
for nonzero v & Ai such that either 



\v\\L2{n) 



(3.1) 



or 



C'{B{v)) < Ce^A'^L 



(3.2) 



holds. Here, C is a constant independent of e and A. The same statements 
hold with Ei and Ai in place of Ei and Ai respectively for i = 2,3. 

The spirit of the proof is similar to that of an unpublished result of Ball 
|Ba06j with the difference that here we take into account the size of the set B 
where \uy\ > 1. According to the theorem, n = is a strict local minimizer 
in the sense that the strict inequality Ei{u) > Ei(0) holds whenever ||u||^2(f^) 
or C'^{B{u)) is sufficiently small and provided n ^ 0. However, as we shall 
see in Section HI n = is a degenerate local minimizer of the Ei in the sense 
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that sets B of positive measure can be introduced at arbitrarily small cost 
measured in terms of Ei{v) — Ei{0) . It would be interesting to determine 
whether, in the terminology of Ball et al [BKM] , [BM84j . u = lies in an 
energy well of the Ei with respect either to ||v||2 or C'^{B{v)). 
We recall from ()2.8p that for any v ^ A 

I ""^^ '^^) ' (/ '^^) ' - ^ ^^'^^ 

for a.e. x G [0, L] where the obvious lower bound dy > C^{lx H B) has 
been used. 

We wish to prove Ei{y) — -E'i(O) > provided v ^ Q and either dx 
or C?(B{uf) is small enough. It is clear that if either 



l^v^yd^y^C^B) 



(3.4) 



or 



vldx>A£.^{B) (3.5) 

n 



and 



then there is nothing to prove. There is also nothing to prove should C^{B) = 
0, since in this case the assumption w 7^ implies (iin view of the boundary 
conditions) that e^^yy + dx > 0, and hence that Ei{v) — -Ei(O) > 0. 
Thus we can assume that 

l^vlydx<^C\B) (3.6) 

/ vldx< AC'^iB), (3.7) 
Jn 

where C^iB) > 0. 

The claim will be that the two conditions (|3.6p and (j3.7p imply lower 
bounds on f^v^dx and C'^{B{u)), thereby proving the contrapositive of 
Theorem 13.11 The intuition behind the claim is explained in the course of 
the next few paragraphs, which should be regarded as a preparation for the 
proofs of Lemma 13.11 and Proposition 13.11 It is on these results that the 
proof of Theorem 13.11 ultimately rests. 
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3.1 Preparation for proof of Theorem 13.11 

Suppose V ^ satisfies (j3.6p and (|3.7|) . Applying Lemma [XT] to (|3.6p . and 
bearing in mind that C^{B) > 0, we see that 

< (3.8) 



t{l-t) 

where 



Note that this automatically implies A > 16e^, which is not a restriction 
since we already know that u = is a global minimizer of Ei in Ai when 
A < 16e^ by Proposition 12.21 Rearranging (j3.8p gives 

A - CHU{B)) - ^" A ' ^^-^^ 

where, as usual, B = B{v) for short. The upper bound merely improves 
slightly the trivial inequality t < 1; the lower bound is new information on 
the set B{v) and is a direct consequence of the assumption (j3.6p . Inequality 
(j3.9p can be interpreted by supposing for the sake of argument that it ap- 
plies to some rectangle B whose sides are parallel to the coordinate axes. Its 
'height' would then be bounded below by a fixed constant. For more gen- 
eral sets B condition eqrefproportionall should be interpreted in an average 
sense, viz. 

The following lemma establishes an inequality involving all the terms 
appearing in the energy Ei{v) and the L^— norm of v. It should be regarded 
as the backbone of Theorem 13. 1[ 



Lemma 3.1. Let v € , let M > and define 



UMiB) = |x G U{B) : ^ v^yyix, y) dy < m| 
Bm = {{x,y)(£B: x G Ha^^)}. 



(3.10) 



Then there is a constant C > independent of v and the dimensions of Q 
such that 



In J \Jn 



t;^(x,y)dx / V(a:,y)dx > U— -^^^^ (3.11) 



M \i:\JiM{B)) 
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Proof. Let x' € IIm{B) and < y < 1. By applying standard results from 
the theory of Sobolev functions we may assume without loss of generality 
that X t—i- is weakly differentiable, and hence that 

v{x, y)va:{x, y) dx = ^v^{x', y) 

on using the boundary condition u(0, y) = 0. The left-hand side of this 
inequality is trivially bounded above by 

C/(x',y):=y^ v\x,y)dx + ^ vj{x,y)dx (3.12) 

for all non-zero a; the dependence on a will be minimized out later. By 
(13.31) we have 



Jl, Ji , vyy dy 



which in view of the fact that x' € IIm{B) implies 



/ v^dy>^C\l,,nBf 



Integrating both sides of this expression over I\.m{B) and applying Jensen's 
inequality to the right-hand side yields 



/ / -'dydx'>^^^^^ (3.13) 

The expression on the right is almost the desired lower bound; a factor of 
C^{IlMiB)) is missing from the denominator. But we know that 

[ [ U{x',y)dydx' > [ [ v^dydx', (3.14) 
Jum{b) Jo JuuiB) Ji^, 

and so it remains to estimate the left-hand side of this inequality from above 
to see that the missing factor can be recovered. Now 



/ / / v'^ {x , y) dx dx' < { v^{x^y)dydx \ dx 

JriMiB) Jo Jo JTiMiB) \Jo Jo J 

= C\I1m{B)) [ v^dx 
Jn 
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The term in can be estimated in the same way. It fohows that 

/ [' U{x',y)dydx'. 
JllMiB) Jo 

Putting this together with (j3.13p and (j3.14p gives 

which on taking = (/^j v^. dx) 2 ( v"^ dx)~ 2 concludes the proof. □ 

We continue to suppose that v G satisfies (13. 6p and (j3.7p . Now by 
definition of 11^^(5) it is the case that 



Mxu{b)\iIm(b){x) < j vlydy, 
from which it follows by integrating and then applying (|3.7p that 

C\U{B) \ Um{B)) < (3.15) 

Hence 

C^{B\Bm) = [ C^{l^r\B)dx 
JnB\UM{B) 

< C\U{B)\Um(B)) 

< • (3-16) 
We are free to choose M = 2Ae^^, thereby ensuring 

C\Bm) > lc\B). 

Now we combine these observations, an upper bound on '^^^^^ and 
Lemma IS.ll to give lower bounds on u^dx and C'^{B). 
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Proposition 3.1. Let v G Ai satisfy (j3.6p and ()3.7p . Then there is a 
constant C independent of v and the dimensions of Q such that 

> \ (3.17) 

(£2(5))2A2 

Furthermore, provided is sufficiently large, 



v^dxj' > (3.18) 

JO-\B{v)) > Ce^A-^L~\ (3.19) 

Proof. Applying Lemma l3.ll to v with the choice of M made above and by 
using inequalities (j3.7p and (j3.9p we see that 



^ 1 C / C^iB) 



(J^„^.2)i Ae-2 yC^iUiB))) 

> 



Ce^ /e^V 1 



1 . I 

l2 



(£2(5))2A^ 

This inequality is (I3.17p . The constant C changes from line to line but it 
remains independent of e, A and L. 

To prove (I3.18P we let p := ||w||l2(q) and q := C^{B)2 and note that 
p.l7p implies 

pq>Ce^A~i. (3.20) 

For brevity we denote the right-hand side of this inequality by /. 

Next, we use the simple interpolation inequality (j2.10p together with 
(13:71) to get 

Ce^A > 

In 

Hence 

p > gq, (3.21) 

where g = CeA~2 . 

Finally, (j2TT]) and ([SZD together imply 
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which in terms of p and q can be written 

p < CLA^q. (3.22) 

The aim is to determine the {p, q) region which is compatible with these 
inequahties. This can be done by looking at Fig [2] below. Note that the line 
with equation p = CLA^q lies above the line with equation p = gq provided 
A > which, in view of Theorem 13.11 is true whenever is sufficiently 
large. This ensures that the required {p, q) region is nonempty. 

It is immediate that solving pq = f and p = gq yields the smallest 
possible value Pmin of p consistent with (j3.6p and (|3.7p . The result is 

Pmin = Ce^A"'^, 

giving the lower bound on ||w||L2(r2) stated in (j3.18p . Similarly, the smallest 
value gmin of q consistent with (13. 6p and (13. 7p is found by solving for q in 
pq = f and p = CLA^q. The result is 

^min = Ce^A-'^L-^ 

giving the claimed lower bound on C'^{B). This concludes the proof of 
Proposition 13.11 □ 

We now draw the preceding results together. 

Proof of Theorem \ 3.1\ We prove the contrapositive under the assumption 
V ^ 0. Suppose Ei{v) < -Ei(O). Then in particular both 

and 

/ vld^<AC\B) 
Jn 

must hold. These equations are (|3.6p and (|3.7p respectively, where, in view 
of u 7^ 0, it can be assumed that C^i^B) > 0. It follows from Proposition ETT] 
that both 11^11^2(^2) > r(e,A) and C^{B) > s{e,A) must hold, concluding 
the proof of Theorem 13.11 

Remark 3.2. The lines p = gq and p = CLA^q coincide when A ~ Ai 
and their relative positions are reversed provided ^ is small enough. Under 
these circumstances the region compatible with all three inequalities (j3.20p . 
p.2ip and (j3.22p becomes empty and the starting point for these calcula- 
tions, namely the inequality Ei{v) < £'i(0), is contradicted. But this makes 
sense since for small it is the case that Ei{v) > -E'i(O) for all non-zero v. 



Local minimizers and low energy paths 23 




9min 



Figure 2: The region compatible with ()3.6p and ()3.7p lies within the wedge 
formed by the lines aibi and 0262 and above the curve containing the points 
ai and 02- 

Remark 3.3. If the lower bounds on ||v||/^2(q) and C^{B) were optimal in a 
scaling sense then as A ^ Ai we might expect r(e, A) — > 00 and s(e, A) — > 
>C^(Q). In other words, the constraints on ||'u||2 and C'^{B{v)) should become 
redundant as A approaches Ai, since when A ^ Ai the state u = is the 
unique global minimizer and the desired inequality Ei{v) > Ei{0) should 
hold for any admissible v. But it is easily checked that A ^ Ai implies 
neither r(e,A) 00 nor s (e,A) C^{n). Thus r(e,A) and s(e,A) would 
appear to be smaller than they could be, i.e. suboptimal. Said differently, 
the inequalities ()3.ip and ()3.2p are likely to be sufficient but not necessary 
conditions for the inequality Ei{v) > -Ei(O). 

4 The effect of modifying the surface energy term 

Assuming for argument's sake that is large enough to ensure, with the 
aid of Theorem 12. 11 that u = is not the global minimizer of Ei for i = 1,2,3 
then it is clear that, regardless of the dynamics governing the passage from 
It = to the global minimizer, the appearance of sets B of positive measure 
is inevitable. Once such sets have been introduced it ought to be possible 
to lower the energy by enlarging them, or by allowing them to evolve in 
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some way in order that the energy is driven to its global infimum. We wish 
to understand the mechanism behind the introduction of sets B of positive 
measure into Q using elements of Ai for i = 1,2, 3. We suppose that such a 
mechanism favours small energy. 

More specifically, we are interested in the least value of Ei{v) — Ei{0) 
consistent with the appearance of sets B of positive measure. It is clear 
that unless one restricts the value of C^{B) the infimum of the difference 
Ei{v) — Ei{0) over Ai will merely reflect the energy of the global minimizer of 
Ei in Ai- The results of Section [3] on local minimality imply that if C^{B) is 
sufficiently small than Ei{v) > Ei{0). In particular, therefore, the quantity 

liminf inf{£;i(t;) : v e Af such that C^{B{v)) < /i} (4.1) 

should be nonnegative, and if it were strictly positive then we could interpret 
it as a lower bound on the cost of 'nucleation'. This is investigated below. 
It turns out that these costs are in fact zero, even in the case of E^ where 
the surface energy includes the full second gradient of v. According to the 
remarks in the introduction, this may not be because the potential W is 
not strongly convex in a neighbourhood of the origin. The force of these 
remarks was that the potential W can be bounded below by a strongly 
convex function in a neighbourhood of zero via a simple argument involving 
a Poincare inequality. Instead, it could be that the only way to ensure 
that the quantity appearing in (14. ip is strictly positive is to replace Wa 
with a strongly convex potential whose second gradient is bounded below 
independently of e and A. 

It happens that low energy competitors can be constructed directly in 
the case of the functional Ei: see Proposition 14.11 below. The construction 
of low energy competitors in the case of E2 and E3 is indirect and is given 
in Section [42] below. 

4.1 Low energy paths for the functional Ei 

Proposition 4.1. Let E'^ be as above and assume is large enough to 
ensure that u = is not the global minimizer of Ei in Ai- Then 

liminf infill (v) : v e Af such that C^{B{v)) < /x} = 0. (4.2) 

In other words, sets B(v) of positive measure can be introduced at arbitrarily 
small energies, as measured by Ei. 
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Proof. The proof of (|4.2p is by construction as follows. Let a, 6 > 0, X > 1, 
f{x) = X and define v:[L-5,L\x [0, 2a] ^ M by 



v{x,y) 



a L - 6 < X < L, 0<y <a 
af{x) - /(^)(^"-^^)' if L - (5 < X < L, a<y<2a. 



Extend v to [L — 5,L] x [0,4a] by reflection in y = 2a, i.e. v{x,y) := 
f (x, 4a— y) if 2a < y < 4a. Finally extend v by zero outside [L—6, L]x[0, 4a]. 
It can then be checked that 

\ a d 

where C is a constant independent of a, 5, e, A and L. It is clear that 
C'^{B{v)) > for each positive 5 and a and for each A > 1. In fact 

C^{B{v)) = Aa5{l - X-^^f , 

so given > we can, by choosing A — 1 sufficiently small and positive, 
ensure that /i > C'^{B{v)) > independently of the choice of a and 5. 
Taking a = 5^ we see that 

Ei{v) < 2CX^5^. 

The conclusion of the proposition follows by letting 5^0. □ 
4.2 Low energy paths for the fuU second gradient functional 

In this section it will be convenient to work with the functional E^. The 
same results then hold for E2 by using the relation (j2.ip . We seek a sequence 
u^-'^ of functions in ^3 satisfying 

(i) > for all j, 

(ii) C^{B{v^j^)) ^ as j ^ 00, and 

(iii) ^3(7;(J)) - ^3(0) ^ as j —>■ 00. 

One way to do this is to take advantage of the fact that if / G and 
if z is the Newtonian potential of / (see e.g. |GTt Chapter 4]), then ZXz = / 
(where A denotes the Laplacian operator, as usual) and 

/ \D^z\'^dx= [ l/pdx. (4.3) 
Jn Jq 
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(This is part of the C alder on- Zygmund theorem, see e.g. \GT\ Theorem 
9.9].) When / is sufficiently smooth, for example Lipschitz continuous, the 
representations 



and 



z(x) 



Vz(x) 



1 

2^ 



ln(|x - y|)/(y) dy 



1 x-y 
2tt |x — yp 



/(y)dy 



(4.4) 



(4.5) 



hold. One now has control of both ||D^z||2 and the pointwise behaviour of 
Vz via the function /. The next lemma details an appropriate choice of a 
sequence of functions f^^^ whose corresponding potentials z^^^ may be used 
to satisfy (i), (ii) and (iii) above. In the following we use the convention 
that Br (a) denotes the open ball of radius r in centred on a. 

Lemma 4.1. Let the sequence of L'^{B2{0),M.) functions f^^^ be defined by 



f^'\x,y) 



where 



A 

R 



ifO<R< 2~i 
if 2'i <R<1 
ifl<R<2, 



a 



R 



2 I 2 

X +y ■ 



Let r]{R) be a smooth cut-off function with support in Bi (0) such that r](R) = 

2 

1 ifO < R <1. Let z^^^ be the Newtonian potential of f^^\ Then the constant 
k can be chosen so that 

(a) ||/*--'-*||2 — > as j — > oo, and 

(b) Zy\Q,(S) > (L^ + 1)5 for all sufficiently large j. 

Here, L is a positive constant. Let ip be a smooth cut-off function with 

support in -63(0) and which satisfies tp{R) = 1 if < R < 1. Then the 
2 

functions z^^'^ := ip{R)z^^^ have compact support in -62(0) and they satisfy 



(c) Zy\{),Q) > (L^ + 1)2 for all sufficiently large j, and 
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f'^) Ib2{0) I^^^Pdx ^0 as j ^oo. 
Proof, (a) To see (a) we compute H/^-'^Hl directly. 

(/(^■))2dx 



2-J 2 



27r^^ + 



16 

Now (a) follows easily (and independently of the choice of the constant k). 

(b) To prove (b) first note that each f^^^ is Lipschitz continuous, implying 
in particular that each z^^^ is (7^(52(0), M) and that (g^]) holds with ^2(0) 
in place of 0,. Thus 



27rz(^')(0,0) 



B2(0) 



-|/(^)(x)dx 



2^TTAj / dR-nAj / R'^^+'^ dR - ttAj 
Jo J2-1 



-2'nA. 



irk 
~2' 



Therefore 



y \-^-J - 4 ^J' 
2 



4^-) (0,0) 



which on choosing k = —6{L +1)2, say, and noting that Aj — > as j 
implies ^■'^(0, 0) > (L^ + 1)2 for all j. This is part (b) of the lemma. 

(c) Part (c) follows easily from (b) and the definition of z^^^ given above. 

(d) By noting that 



00, 



|L>2z(j')|dx < C 



B2{0) 



Ua)|2+[Vz(i)|2+|Wj-)|2dx 



B2(0) 



Iz^-')!^ dx and f iVz^-'^l^ 



for some constant C, it suffices to prove that /j 
converge to zero as j 00. The convergence to zero of the term /^^(o) 1-^^-^^''^ I 
is guaranteed by (|4.3p and part (a) above. By the representation (j4.4p . stan- 
dard estimates and Fubini's theorem, 



j^-^') dx < C 



'B2(0) 



(ln(|x-y)|)2dx /(■'■)(y)2dy. (4.6) 



B2(0) I JB2(0) 
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Therefore \ \z^^^\\l as j — > oo. To check the convergence of ||Vz(^')||2 to 
zero, write 

|Vz(^')pdx= / Vz(j') • dzy - / z^Az^^^dx. (4.7) 

Next, note that both z^^^ and Vz*--'^ converge uniformly to zero on dB2{0), 
which can be verified by using the fact that each f^^^ has compact support in 
-63(0) together with the representations (|4.4|) and (j4.5|) . The second term in 

(|4.7p can be estimated by using Holder's inequality and (|4.3p in that order, 
giving 

/ Iz^-'') Az(^')|dx< C||/(^')||^ 

Jb2(0) 

for some generic constant C. Therefore ||Vz^-'''||2 0, which concludes the 
proof. □ 

The next result formalises the statement made at the start of this sub- 
section. The proof can easily be adapted to show that sets B{v) of positive 
measure can be introduced into 0, in such a way that E2{v) — £2(0) can be 
made arbitrarily small. 

Theorem 4.1. Let and A3 be as per ()1.3p and ()1.7p respectively. Then 
there exist sequences {v^^^} C .43 such that 

(i) £'^{B{v^^^)) > for all sufficiently large j, 

(a) C'^{B{v'^^^)) ^ as j ^ 00, and 

(Hi) E3{v^^^) - ^3(0) ^ as j —> 00. 

In other words, sets B(v) of positive measure can be introduced at arbitrarily 
small energies, as measured by E3. 

Proof. Let P be the point |) in 17 and define the planar affine map T by 



Then T{Vt) C ^2(0), hp(T) = ^ and we can define 



^;0■)(x) = ~z^^\t[x)) 
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for X £ 0,, where z^^^ is as per Lemma l4.1i To check that (i) holds is now 
straightforward. Indeed, since 

4^-)(x) = hp(T)4^-)(r(x)) 

for all X it follows from part (c) of Lemma l4.ll that 

vlj>\p) = hp{T)z^^\0) > 2. 

Since u^^) is it follows that ^^({l^-'^l > 1}) > 0, which is statement (i) 
above. To see statements (ii) and (iii) note that 

^3(t,0-)) = [ e2|D\(j')|2 + 4i)^dx + A£2(5(i;(i))) 
Jn 

< [ eHMT)f\D^z(^^\' + m'dy + C [{4^^fdy 
Jb2{o) Jn 

< C I |Z)2zO')|2 + |V5(j')|2dy, 

J 32(0) 

where we have applied Chebychev's inequality in the second line. The right- 
hand side can now be made arbitrarily small by appealing to part (d) of 
Lemma |4.H proving parts (ii) and (iii) of the theorem. □ 

Remark 4.2. It may be significant that none of the examples constructed 
above has small support. It could be that if we require small support (which 
is physically reasonable) then there may well be an energy barrier associated 
with the appearance of sets B(v) of positive measure in 0. It is not clear 
how or why any such condition should be imposed a priori. Indeed, there are 
many possible constraints on such sets which may be physically reasonable 
and yet do not enter into these variational models. Other factors, such as 
the asymmetry of the boundary conditions, may also have had a role to play. 



5 Appendix 

This section is included for completeness only. We show that the global 
minimzer C/j, say, of Ei in Ai exists. Only the case i = 2 is considered here: 
the others follow by analogy. 
Recall that 

E2{v)= [ \V{vy)\^ + vldx + AC\A{v)) 
Jn 
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where 

A{v) = {(x,y) . \yy(^x,y)\ < 1}. 

Now E2 is bounded below by zero, so its infimum in A2 exists. Let {v^^^} C 

A2 be a minimizing sequence. Then, since ^{vy^ ) is bounded in L^, it follows 

that there is w in such that for a subsequence (and after relabeling) 

Vy'^ ^ w in Ty^'^. By the Rellich-Kondrachov compactness theorem we 

(7) 

can suppose that the sequence Vy converges strongly, and is in particular 
bounded, in L^. Looking again at E2 it follows that vji^ is also bounded 
in L^. So iVv*'-'^! is bounded in L^, and hence there is some function U in 
Tyi'2(0,M) such that 

yU) -.u in W^'^. 

It follows that w = Uy. The trace theorems for Sobolev functions now 
imply that U G A2- (One could use [EGj Section 4.3, Theorem 1], or 
(A. 5) in the appendix of [KM94j . for example.) Finally, the sequential lower 
semicontinuity of 

V / f ^ dx 
Jn 

and of 



Jn 



with respect to weak convergence in W'^''^, together with Fatou's lemma, 
imply that 

\iminfE2{v^^^) > E2{U), 
concluding the proof of the existence of the global minimizer of E2 in A2- 
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